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ABSTRACT 



The transverse forces acting on smooth and sand-roughened circular 
cylinders immersed In an harmonically oscillating flow have been mea- 
sured using a recently constructed U-shaped water tunnel. 

The i i ft coefficient and the frequency of the alternating force 
have been determined. For smooth cylinders, the lift coefficients were 
found to depend on both the Reynolds and the Keu legan-Carpenter numbers 
for Reynolds numbers between about 20,000 and 150,000. For rough 
cylinders, the lift coefficients were found to be Independent of the 
Reynolds number within the range of relative roughnesses and Reynolds 
numbers tested. The results have also shown that the lift coefficients 
for smooth cylinders at low Reynolds numbers are nearly identical with 
those obtained for rough cylinders at very high Reynolds numbers at 
the corresponding amplitude ratios. For both rough and smooth cylinders, 
the transverse force is a significant fraction of the in-line force 
and must be taken into consideration in the design of structures. 

it is recommended that the experiments be extended to the self- 
excited hydro-elastic oscillations of cylinders in harmonic flows. 



J00H32 iU'j-HA. 



TABLE OF CONTENTS 



I. INTRODUCTION 10 

II. ANALYSIS OF THE DATA AND THE FORCE COEFFICIENTS 14 

A. GOVERNING PARAMETERS . 15 

III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 19 

A. U-SHAPED OSCILLATING- FLOW TUNNEL 19 

B. CIRCULAR CYLINDER MODELS 25 

C. FORCE MEASUREMENTS 26 

0. ACCELERATION, ELEVATION, OR VELOCITY MEASUREMENTS 31 

E. DATA REDUCTION 35 

IV. DISCUSSION OF RESULTS 36 

A. BLOCKAGE AND LENGTH-TO-D I AMETER EFFECTS 36 

B. TRANSVERSE FORCE DATA 38 

1. Smooth Cylinders — 38 

2. Sand-Roughened Cylinders R 2 

V. CONCLUSIONS AND RECOMMENDATIONS 65 

APPENDIX A Computer Programs- 67 

APPENDIX B Transverse Force Data for Smooth Cylinders — - 69 

APPENDIX C Transverse Force Data for Rough Cylinders — 82 

LIST OF REFERENCES 105 

INITIAL DISTRIBUTION LIST 106 



) 



i 

1 

a 

i 

i 

i 

i 

) 

\ 

i 

i 

i 

J 

h 

) 

* 

t 

i 

! 

4 

3 

; 

| 



} 



4 

1 

l 



5 






LIST OF FIGURES 



Figure 

1. Schematic drawing of the U-Tunnel- ------------ 20 

2. Photograph of the U-Tunnel ---------------- 21 

3. Photograph of butterfly valve system 23 

4. Scanning Electron Microscope photographs of 

sand-roughened surface 27 

5. Gage assembly- ---------------------- 29 

6. Acceleration and transverse force traces --------- 32 

7. Position of pressure taps- - — ------------- 33 

8. C^(max) versus K for 8 3 497 - — - ------- — 39 

9. C^(max) versus K for 6 3 784 --------------- 40 

10. C^(max) versus K for 6 - 1107- 41 

11. C^(max) versus K for 8 = 1985- 42 

12. C^(max) versus K for 8 3 3123 - - — - — ------ 43 

13. C^(max) vr.'sus K for 8 3 4480- 44 

14. C^(max) versus K for 8 ’• 5260 45 

15. C^(max) versus K for various values of K and 8 ----- - 47 

16. C^(max) versus Re for various values of K- ------- - 48 

17. Relative frequency of vortex shedding as a function 

of the Reynolds and Keu legan-Carpenter numbers - - - 50 

18-21. C. (max) versus K for k/D = 1/370 and 

various values of 8- ----- — 53 

22-25. C (max) versus K for k/D = 1/195 and 

various values of 8- ------------------- 57 

26. Combined plot of C. (max) data for k/D = 1/195 

for various va I ues of 8- - - - 61 



6 



NOMENCLATURE 



A Virtual amplitude of the motion 

Aj Amplitude of the motion 

Maximum acceleration 
m 

CLMAX Maximum transverse force coefficient, also usea as C^(max) 
CLRMS Root-mean-square transverse force coefficient 

0 Diameter of the test cylinder 

F instantaneous total force acting on the test cylinder 

f Frequency ratio, f y /l/T 

f Frequency of the first harmonic of the transverse force 

F Measured force 

ni 

g Gravitational acceleration 

H Distance between pressure taps and mean water level 

K Keuiegan-Carpenter number 

L Length of the test cylinder 

Re Reynolds number 

k Sand roughness particle size 

s Distance between the pressure taps 

T Period of oscillations 

t Time 

U, Maximum velocity 

• m 

u Instantaneous velocity 

w Width of the test section 

8 Frequency parameter 

Y Specific weight of fluid 

♦ 
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Ap 

sr 

V 

P 



Di f f erentl a I pressure 
Phase angle 

Fluid kinematic viscosity 
Fluid density 
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I. INTRODUCTION 



Information about the time-dependent forces acting on bluff bodies 
In general and circular cylinders In particular has considerable practi- 
cal interest In ocean and wind engineering, as well as In basic under- 
standing of fluid mechanics. Extensive discussion about the flow-induced 
forces and oscillations exists in the literature, but the basic hydro- 
dynamic data are lacking, particularly at higher Reynolds numbers which 
are of current practical Interest* 

Much ot the present knowledge on the hydrodynamics of the oscillatory 
flow about the cylinders has been obtalnud by means of model tests In 
wavo channels at Reynolds ni:r,iDers generally two to three orders of magnl- 
iude smaller than prototype Reynolds numbers* These model tests, which 
have been primarily concerned with in-line forces, have disclosed that 
the drag and Inertia coefficients may be correlated with the amplitude of 
the motion relative to a characteristic dimension of the oody and that 
the effect of the Reynolds number is obscured by that correlation. The 
field data have shown rather large scatter and no systematic correlation 
with either the relative amplitude or the Reynolds number. Understand- 
ably, there Is considerable controversy In the literature regarding the 
effect of viscosity, free surface, etc. 

Vortex shedding gives rise to an alternating force In both steady and 
time dependent flows. In spite of the considerable Interest, however, 
the transverse force or the lift force In harmonic flow received very 
little attention. Recently, It became clear from the observations of the 
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oscillations of long piles and strumming of cables thct the lift forces 
are Important net only because of their magnitude buT also because of 
their alternating nature which under certain circumstances may lead to 
the phenomenon known as lock-in or vortex synchronization. This phenomenon 
may cause failure due to fatigue and Increased in-line force. Obviously, 
the total instantaneous force acting on the structure Is Increased by 
the lift force and modified by the oscillations of the body. This In- 
crease refers to the vectorial sum of the in-line and transverse forces 
and not to the forementioned oscillations In the In-line force due to 
vortex shedding. 

Some of the previous studies Include those carried ou+ by Chang Cl], 
Bidde C2D, Wiegei and Delmonte C3], Mercler Sarpkaya and Tuter C5], 

Jones Q6H, Isaacson [7], and Sarpkaya C8]. Bidde [ 2 ] dealt primarily 
with the ratio of the transverse force to the In-line force in many flows 
and concluded that the lift force behavior is primarily dependent on the 
Keulegan-Carpenter number rather than upon the Reynolds number and that 
the predominant lift frequency Is twice the wave frequency. Bldde f s 
data are difficult to Interpret because of the fact that the Reynolds 
number and the Keulegan-Carpenter number were calculated In terms of some 
average value, that the force measured was the total force on the complete 
length of the pile, and that the submerged end of the vertical cylinder 
was completely free to generate a complex three dimensional flow and 
Influence the vortex shedding. 

Wiegei and Delmonte C3] extended Bidde f s work and used the Keulegan- 
Carpenter number based on the wave-surface kinematics. They have in 
general confirmed Bldde’s conclusion except for the fact that the lift 
frequency was irregular and varied from about 1.3 to 6 times the wave 
frequency. 
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Sarpkaya and Tu+er C5D measured the transverse force on cylinders In 
plain harmonic flow at relatively low Reynolds numbers and found that 
the maximum lift coefficient Is primarily a function of K and that It 
can acquire large magnitudes near Keulegan-Carpenter numbers of 15. This 
work was subsequently extended to cylinders In the vicinity of a plain 
wall C63 and to higher Reynolds number flows [8]. 

Isaacson measured the lift force on vertical cylinders In wavy 
flows within a Keu legan -Carpenter range of about 0 to 25 for Intermediate 
depth waves. The Reynolds number range covered was from about 100 to 
5000. Isaacson concluded on the basis of his and others work that lift 
was dependent both on the Reynolds and Keulegan-Carpenter numbers and 
’hat the dependence of the lift coefficient on the Keulegan-Carpenter 
number was considerably stronger and tends to obscure the weaker depend- 
ence on the Reynolds number. He also argued that for higher values of 
the Keulegan-Carpenter number the predominant lift frequency must In- 
crease with the Keu legan-Carpenter number. 

» 

It Is apparent from the foregoing that both the In-line and transverse 
forces on the structures in the marine environment must be evaluated with 
proper attention to the order of magnitude of the Reynolds and Keulegan- 
Carpenter number**. Furthermore, one must take Into consideration the 
surface conditions of the structure such as rigid and flexible marine 
growth. There Is not at present either sufficient or reliable data to 
answer the needs of the designer and to isolate the influence of the Indi- 
vidual factors on the variation of the force transfer coefficients. In 
view of the foregoing considerations, the present research program was 
undertaken with the following objectives: 



(a) To furnish data, obtained under carefully controlled laboratory 
conditions, about the in-line and transverse forces acting on 
smooth cylinders in an harmonically oscillating fluid at Reynolds 
numbers up to 1,000.000 and 

(b) To identify the physical mechanism and parameters responsible 
for the correlation or scatter of the force-transfer coeffi- 
cients. 

This thesis deais only with the transverse forces acting on smooth 

t 

and sand-roughened circular cylinders, it does not deal with ocean waves, 

\> 

non-harmonic fluid oscillations, wave and current combinations and its 

» 

consequences, diffraction effects, free-surface and/or wall proximity 
effects, fluid elasticity or hydro-elasticity of flexible or fiexibly 
supported cylinders in harmonic motion, or finally the interference ef- 
fects between neighboring structural elements. 

The transverse-force coefficient for both smooth and rough cylinders 
have been obtained experimentally and expressed in terms of Reynolds 
number, Keu iegan-C~rpenter number, and relative roughness. In addition, 
the relative frequency of vortex shedding has been determined as a 

4 

function of the same parameters cited above. 
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I I . ANALYSIS OF THE DATA AND THE FORCE COEFFICIENTS 



A completely satisfactory analysis of the resistance In unsteady 
separated flow has escaped the concentrated efforts cf many researchers. 
No theoretical model can, at present, predict the complete force and 
flow character! sti cs of an harmonic motion about a circular cylinder. 

In the absence of such an analysis, the most serious difficulty lies in 
the description and i nterpretat ion of the history of the flow and the 
effect of vortices and turbulence in the flow approaching the body at a 
particular time in the oscl Nation. 

in the present study the following transverse force coefficients 
have been defined and evaluated: 

CLMAX “ ^ ay i mum P eak the transverse force in a cycle ^ I j 

(0.5)(p)(D)(L)(U m 2 ) 

and 



Clrms - root-mean-square value of the transverse force in a cycle 

O.SpDLU 2 
m 



( 2 ) 



The frequencies of oscillation have been expressed either in terms of 
the frequency of oscillation of the harmonic flow as 



_ frequency of the f I rst harmonic of the lift force in a cycle 
r frequency oi 1 the harmonic oscillation 

(3) 



or in terms of the Strouhal number defined by 



St = 



f v D 




f 

r 

K 



(4) 



i 

i 



in which p represents the density of fluid, 0 and L the diameter and 
length of the cylinder, U m the amplitude of the oscillating velocity and 
f v the frequency of the transverse forces. 

A. GOVERNING PARAMETERS 

The coefficients cited above v/iil have to be correlated through the 
use of suitable parameters in order to show that they have some degree of 
universality within the range of the parameters encountered. A simple 
dimensional analysis of the flow under consideration (uniform harmonic 
motion about a rough circular cylinder placed with Its axis normal to 
;he flow) show that the time-dependent force coefficients may be written 
as: 



f( 



U T 
m 



O.SpDLU 



U D . . 

m J< t. 

v ' 0' T 



■ f(K, Re, £, i> 



(5) 



(5a) 



Evidently, U m T/D Ma y replaced by 2nA/D or simply by A/0 where A re- 
presents the amplitude of the oscillations. 

There is no simple way to deal with equation (5a) even for The most 
manageable time-dependent flows. The evaluation of the instantaneous 
values of the force coefficients is not at present feasible. Another 
and perhaps the only other alternative is to eliminate time as an inde- 
pendent variable in equation (5a) and consider suitable time-invariant 
averages or amplitudes of the force coefficients. Thus, one has 



CLMAX 

CLRMS 



U T . 

, m _ k. 

( n > 



(6) 



St 






Equation (6), in spite of its simplicity, gives rise to many ques- 
tions: Do the averaged force coefficients really depend upon both Re 

and K; Are K and Re the most suitable governing parameters; Can one ob- 
tain meaningful conclusions by plotting the data for a given coefficient 
with respect to, say, K and connecting points having equal Re or vice 
versa; How should the experiments be conducted so that equation (6) 
yields manageable plots; Which of the two parameters has a more pro- 
nounced effect on the force coefficient under cons i derat Ic:i; Why has 
there been considerable scatter in the field data when plotted with 
respect to either K or Re; Are there ranges of K and Re in which the 
effect of one is obscured by a reasonable correlation of the force co- 
efficients with the other? These a~d similar questions have been raised 
by many other investigators and attempts were made to establish suitable 
correlations. The state of the art is such that the past conclusions 
and conjectures can be critically scrutinized only through the acquisi- 
tion of data of high intrinsic quality obtained under controlled labora- 
tory conditions with relatively simple harmonic flow situations. The 
purpose of such an effort i s by no means to remove the need for actual 
full scale experience, in fact it is to encourage full scaie experiments 
and to enable those concerned to interpret and better understand the 
factors affecting the force-transf er function. 

Past experience [9, 5} has shewn that the force coefficients obtained 
under controlled laboratory conditions are primarily functions of K at 
relatively small Reynolds numbers (Re < 25,000) and that the effect of 
viscosity is obscured by the fairly good correlation between K and the 
force coefficients. Again previous efforts and the reasoning based on 
dimensional analysis have shown that there is an undeniable effect of 



\ 

A 



the Reynolds number. Thus, means have to be devised to delineate the 

effect of both K and Re or some other viscosity dependent parameter [8]. 

It Is apparent that the maximum velocity appears In both 

K = U T/D and Re = U D/v. Simple rules of dimensional analysis state 
m m 

that one obtains the maximum amount of experimental control over the 
dimensionless variables, if the original variables, that can be regulated, 
each occur in only one dimensionless product. Thus, If U m Is easily 
varied experimentally, then U m should occur In only one of the independ- 
ent dimensionless parameters. With this hint in mind. Re in equation (6) 
Is replaced by Re/K = 0 /vT. This parameter shall be called the 
-frequency parameter- and denoted by 0, so that 

6 * D 2 /vT * (7) 

Evidently, for a series of experiments conducted with a cylinder of a 
given diameter in water (of uniform and constant temperature and density) 
undergoing harmonic oscillations with a constant period (T), 0 Is held 
constant. Then the variation of a force coefficient with K may be 
plotted for constant values of 0. Subsequently, one can easily recover 
the Reynolds number from 



Re = 0K (8) 

and connect the points, on each 0 a constant curve, representing a given 
Reynolds number for a family of suitably selected values of the Reynolds 
number. Such a procedure eliminates the difficulty of trying to draw 
contours of constant K, or constant Re, or constant CLMAX or CLRMS versus 
K or Re or K versus Re. Suffice It to note that the data reported here-ln 



shall be analyzed according to the relationship 



C, = f,(K, 6, £> (9) 

Ca coefficient} 

and the Reynolds number will be used in the manner cited above. The 
power of this new plotting procedure (new, as far as the wave force 
analysis is concerned) will become apparent later. 



III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 



A. U-SHAPED OSCILLATING- FLOW TUNNEL 

Experiments carried out by Sarpkaya and Tuter C5j in the initial 
phases of the study with small smooth cylinders at relatively low Rey- 
nolds numbers have proved the versatility and the usefulness of a U-shaped 
osci I lating-f low apparatus. Thus, in an attempt to achieve larger Rey- 
nolds numbers, it was only natural to construct a larger U-shaped tunnel. 

Among the various designs, the one shown in figure I was finally 
selected for construction [8], A photograph of the completed and fully 
instrumented structure is shown in ftgure 2. It conststs of eleven 
modules for ease of construction, transportat Ion, and mounting. Each 
module is made of 0.95 cm. (5/8 inch) aluminum plates and reinforced with 
1.27 x 10 x 46 cm. (1/2 x 4 x 18 inches) aluminum flanges welded to the 
plates. The modules were assembled with the help of an air drying sili- 
con rubber seal between the flanges of two adjacent modules and 2.54 cm. 

(I Inch) steel bolts placed 15 cm. (6 inches) apart. The inside of each 
module was precision machined so that the largest misalignment was about 
I mm. (0.04 inches). 

Prior to the description of its instrumentation and operation, a 
few words are necessary aboul the general shape of the tunnel. The 
cross-section of the two legs is 183 x 91.5 cm. (6x3 feet). This 
selection was dictated by several considerations such as available ceil- 
ing height, pressures to be encountered and hence the structural and 
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Figure 2. Photograph of the U-Tunne/ 
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economic considerations, desire to obtain a virtual amplitude or 
velocity of oscillation at least twice that of the free surface, period 
of oscillation, Reynolds number and the relative amplitude A/D desired, 
natural -damping of the oscillations, and the magnitude and the fre- 
quency of the forces. The length of the horizontal test section was 
chosen larger than twice the virtual amplitude to insure fully developed 
uniform flow at the test section. Finally the two corners of the tunnel 
were carefully streamlined to prevent flow separation (see figure* I). 
This design proved to be more than adequate for no separation was en- 
countered, and also the desired frequency and amplitude of oscillation 
were achieved. 

The auxiliary components of the tunne* consisted of plumbing for 
the filling and emptying of the tunnel with hot and cold water, a heat 
exchanger, butterfly-valve system, and the air supply system. The 
plumbing consisted of simple piping, valves, a small pump, and a filter. 
None of these will be described further. 

The butterfly-valve system (mounted on top of one of the’ legs of 
the tunnel) consisted of four plates, each 4.6 cm. (13 inches) wide 
and 92 cm. (36 Inches) long. A 2.54 cm. (I Inch) steel shaft was placed 
at the axis of each valve plate. (see figure 3). Aluminum housings 
Supported both ends of the shaft with sel f-al igning ball bearings, A 
15 cm. (6 inches) gear was attached to one end of each shaft which ex- 
tended beyond the bearing. All four valve plates were then aligned and 



# 

The vii „al amplitude is defined as the amplitude of oscillation 
which the cylinder experiences. Here it is exactly twice the amplitude 
of oscillation of the free surface in one leg of the tunnel. 
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driven by a simple rack and pinion system. The rack was actuated by an 
air-driven piston with the help of a three-way valve connected to the 
laboratory air supply system. 

The valves, in their closed position, completely sealed the top of 
one of the legs of the tunnel. The top of the other leg was left open. 
Initially, the butterfly valves were closed and air was admitted to that 
side of the tunnel to create the desired differential water level between 
the two legs of the tunnel. Then the valves were opened quickly with the 
help of a pneumatically-driven three-way control valve. This action set 
the fluid in the tunnel in oscillatory motion with a natural period of 
T - 5*507 seconds. A series of experiments was conducted with one of 
the test cylinders to evaluate some of the experimental characteristics 
of the tunnel. It was found that the damping of the motion is such that 
the amplitude of osci I lations .decreases about 3 mm. (0.12 inches) per 
cycle and about 0.15 mm. (0.06 inches) per cycle for amplitudes half or 
less than of the maximum. Thus, over a period of four or five complete 
cycles of oscillation at any amplitude, the amplitude, velocity, and 
the acceleration of the fluid changed about one percent. Evidently ^he 
forced oscillations of the fluid; if such a methoo were 1*0 be employed, 
cannot yield the amplitude to an accuracy better than one percent. 
Additionally in such a method one has to contend with some nigh frequency 
vibrations, however smail, superimposed on the acceleration. These re- 
sult from the cyclic operation of the butterfly valves. It is because 
of these considerations that the experiments were carried out by letting 
the system damp out the amplitude over many cycles of oscillations. The 
advantages of the method adopted become apparent very quickly. Firstly, 
the oscillations were so smooth and quiet that one could not know or 
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even hear that 5000 gallons of water were in oscillation. The elevation, 
acceleration, and all force traces were absolutely free from secondary 
oscillations so that no filters whatsoever were used between the trans- 
ducer output and the recording equipment. Secondly, the method adopted 
enables one to cover ail possible values of K for a given B and k/D and 
see the evolution of the forces over a period of about 30 minutes. 

3. CIRCULAR CYLINDER MODELS 

Circular cylinders with diameters ranging In size from 20.3 cm. (8 
Inches) to 5 cm. (2 Inches) have been usad in this study. The cylinders 
were turned on a lathe from aluminum pipes or plexiglass rods and 
polished to a mlrror-shlne surface. There is no doubt that the surfaces 
were hydrodynamical ly smooth. The length of each cylinder was such that 
it allowed 0.08 cm. (1/32 Inch) gap between the tunnel wall and each end 
of the cylinder. As will be noted later, the cylinder was prevented 
from moving towards one or the other wall by means of small 0-rings 
attached to the round cantilever end of the force transducers. A double- 
ball precision bearing (SKF-2303-J) with an approximately 1.5 cm. (0.6 
Inches) bore was Inserted at each end of the cylinder in aluminum hous- 
ings which sealed the cylinder air tight. The other face of each bear- 
ing was flush with the end of the cylinder. The same cylinders were also 
used as rough cylinders by coating them with sand of desired size. For 
this purpose clean sand was screened through +he use of the combination 
of appropriate size sieves to separate the desired size. The cylinder 
surface was then cleaned with alcohol and carefully coated with an epoxy 
paint. After considerable practice of application of the epoxy paint 
it has been possible to obtain an extremely thin and uniform coating 
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which was free of ridges, waves and local buildup. Next, the selected 
sand was sprinkled over the gently- rotating cylinder with a sieve 
slightly larger than the sand particles. The excess sand was easily 
removed after about 15 minutes through gentle rubbing by hand. The re- 
sulting roughness was perfect beyond expectation as evidenced by the 
two sample photographs shown In figure 4. The pictures of the rough 
surface were taken with a scanning electron miscroscope. The actual 
size of the sample sand Is shown in the figure. 

it was necessary to choose the proper size of the sand for each 

cylinder in order to maintain the k/D ratios (k/D ~ 1/195 and k/D - 

1/370), besides the smooth cylinder, selected* Furthermore, it was 

necessary to carry out the experiments for a given cylinder at various 

2 

temperatures in order to vary D /vT or, in other words, to investigate 
the variation of the lift coefficient with Reynolds number for a given 
Keu legan-Carpenter number and relative roughness. The water in the 
tunnel was heated to the desired temperature through the use of the 
heat exchanger cited previously. The temperature in the tunnel remained 
constant for a period of several hours thus, it was not necessary to 
continuously circulate the water through the heat exchanger. Evidently, 
the enormous heat capacity of the tunnel and the fluid helped to main- 
tain a constant and uniform temperature. The expansion of the test 
cylinders was less than 1/32 inches (0.8 mm) at the highest temperature 
tested. 

C. FORCE MEASUREMENTS 

Two identical force transducers, one at each end of the cylinder, 
were used to measure the instantaneous in-line and transverse forces. 
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Figure 4. Scanning Liectron Microscope Photographs of sand-roughened 

surface 



27 




The basic transducer was manufactured by B.L.H. Electronics, Inc. under 
the trade name LBP-I and catalogue no. 420271. A photograph of the gage 
assembly is shown in figure 5. The gage had a capacity of 2224 N (500 
Ibf.) with an overload capacity of 200 percent. The deflection of the 
gage under 500 Ibf. load was 0.25 mm. (0.01 inch). For the largest 
cylinder and amplitude encountered in the experiments, the maximum load 
was less than 200 Ibf., and the deflection of the beam was less than 
0.2 mm. (0.008 inches). 

A special housing was built for each gage so that it can be mounted 
on the tunnei window and rotated to measure either the in-line or the 
transverse force alone. The bellows which protected the strain gages 
had to be waterproofed in such a manner that they would not adversely 
effect the operation of the gages when subjected to about 6 m. (20 feet) 
water pressure at temperatures i8 C (64 F) to 74 C (165 F). For this 
purpose the bellows were completely filled with Dow Corning 3140— RTV 
coating without bringing the rubber into contact with air during the 
filling operation. After filling, the ends of the bellows were sealed 
ai 1 " tight with special clamps. The silicon rubber remained in its 
original liquid form throughout the operation of the gages. 

The cylinders were placed in the test section by retracting the 
gages from their housing and then pushing them Into the bearings mounted 
on each end of the cylinders. As noted earlier the 0-rings placed on 
the cantilever end of each gage prevented the test cylinders from moving 
sideways towards one or the other wall and helped to exactly set the 
0.8 mm. (i/32 inch) space between the cylinders and the tunnel wall. 

The cylinders w^re free to rotate at the application of a slight torque 
by hand. 



/ 
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Figure 5. Gage Assembly 
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After the mounting of the first cylinder, the exact angular position 
of the gages within their housing had to be determined and set with a 
pin so that the gages measure either only the in-line or the transverse 
force. For this purpose, an approximately 4C0 N load was hung on the 
cylinder with a lubricated nylon rope. The in-line force (acting in the 
horizontal direction) was observed on the amplifier recorder system. 

Then the gage was rotated in small increments until the in-line force 
was exactly zero. A final check was made by measuring +he outputs of 
the gages with a precision voltmeter. The position of each gage was 
marked and set with a pin. Finally, four boits were placed on the gage 
housing to hold the gages rigidly in position. Removal of these bolts 
and the pin allowed the rotation of the gages exactly 90 degrees, after 
which the bolts and pin were placed in position. In this manner the 
gages were capable of measuring either the in-line or the transverse 
force without any "cross talk" between the in-line and transverse forces. 
Ordinarily, one gage was set to measure In-line force and the other gage 
the transverse force. At times both gages were used to measure only 
the in-line or the transverse force. 

The calibration of each gage was accomplished by hanging loads in 
the middle of the cylinder after setting both gages to sense only the 
transverse (here vertical) force. The directional sensitivity of the 
gages was also checked by applying identical loads upwards on the test 
cylinders with the help of a hook-cantilever arm attached to the top of 
the tunnel outside the test section. Repeated calibrations have shown 
that the gages are perfectly linear up to 2000 N; they yield the same 
signal for loads applied either downward or upwards; and that the gages, 
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together with the electronic system to which they were eventually con- 
nected were capable of sensing forces as small as 0.1 N (0.02 I b f . > . 

The in-line and transverse force were simultaneously recorded with 
the instantaneous acceleration on two two-channel Honeywell recorders 
running at a speed of 10 divisions per second. This speed gave 55.07 
divisions per cycle. The amplitude of the transverse force and the flow 

characteristics such as U T/D and Re were determined from these traces. 

m 

RMS value of the lift force wa3 determined for each cycle by reading the 
force at every division or 0.1 second intervals. Sample transverse 
force and elevation traces are shown in figure 6. 



D. ACCELERATION, ELEVATION, OR VELOCITY MEASUREMENTS 

It is because of the extreme Importance of the accurate measurement 
of the instantaneous value of these quantities that they are discussed 
here separately. 

Firsvly, it should be noted that the measurement of the amplitude, 
acceleration, elevation, or the velocity is a matter of i nterpretat ion 
of the signal received from the appropriate transducer in light of one 
of the following expressions 



2ttA _ dU . _ 2 tt 

U m = — a m * at " ( T ) A T U m 



in which T = 5.507 seconds for the experiments reported here. 

Three transducers were used to generate three independent d.c, 
signals, each proportional to the instantaneous value of one of the 
quantities cited above. The first one consisted of a seven foot long 
platinum wire stretched vertically in one leg of the tunnel. The out- 
put of the capacl .ance-wl re bridge was connected to an eight channel 




Figure 6, Acceleration an<l transverse force traces 



amplifier-recorder system. The response of the wire was found to be 
perfectly linear within the range of oscillations encountered. The 
wire was capable of yielding a measurable signal for changes in water 
elevation as small as 0.8 mm. (1/32 inches). Such a sensitivity was 
not, however, always desirable for the instabilities on the water sur- 
face gave rise to small oscillations in the analog records. The 
effect of such insrabl I ities were practically eliminated by placing 
the wire along the axis of 30 cm. < foot) diameter and 2 l 3 cm. (7 
feet) long thin plas+lc pipe. Be that as It may, the use of this wire 
was rendered unnecessary due to the use of a more reliable method. 

The second method consisted of the measurement of the instantaneous 
acceleration by means of a differential-pressure transducer connected 
to two’ pressure taps placed horizontally 61 cm. (2 feet) apart and 
91.5 cm. (3 feet) to one side of the test section. The output of tne 
transducer was again connected to the eight channel recorder. The in- 
stantaneous acceleration was then calculated from ip = ps(dU/dt) where 




Figure 7. Position of pressure taps 
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Ap is the differential pressure, s the distance between the pressure 
taps and dU/dt Is the instantaneous acceleration of the fluid. The 
effect of the pressure drop due to the viscous forces over the distance 
s was calculated to be negligible. 

The third method again consisted of the measurement of the differ- 
ential pressure between two pressure taps. The two taps were placed 
symmetrically on the two vertical legs of the tunnel at an elevation 
127.0 cm. (50 inches), (see figure 7) below the mean water level, i.e., 
H = 38 in. Applying Bernoulli’s equation for unsteady flow between 
each pressure tap and the instantaneous level of water, it is easy to 
show that twice the amplitude of the free surfaco oscillation (virtual 
ampl It de) is given by 



A - 2Aj 



Ap/y 



- - (^) 2 H 
a .T 



(10) 



in which g and T are constant and H is kept constant. Thus the signal 
of this transducer yielded the virtual amplitude cr the maximum 
velocity in each cycle. It was entirely free frvn noise or small free 
surface effect. The transducer was calibrated and Its linearity checked 
before each series of experiments. 

Suffice It to say that al! three methods gave nearly identical 
results and yielded the amplitude, velocity, or acceleration to an 
accuracy of about two percent relative to each other. These comparisons 
as well as the perfectly sinusoidal ard noise-fr9e character of all 
pressure and force traces speak for the suitability of the unique test 
facility used in In is study. 



E. DATA REDUCTION 



Experiments were repeated at least three times for each cylinder. 
Often only two, but at times as many as four runs were evaluated. 
Suitably selected cyles were read every 0.1 second and the tabulated 
data were transferred to the computer data cards together with addi- 
tional relevant information such as calibration factors, cylinder 
diameters, number of data cards, etc. A computer program calculated 
the force transfer coefficients previously discussed (see Appendix A). 

It should be emphasized once again that only the amplitude and 
frequency of the first harmonic of the transverse forces were evaluated. 
Ordinarily one would write the transverse force as 



and attempt to evaluate the coefficients a and b or the coefficients 

n n 

a^ and the phase angle for each harmonic. In view of the fact that 
this Is a highly time consuming affair and that the first harmonic is 
almost always several , times larger than the second or higher harmonics. 
It was decided to evaluate only the character! st I cs of the first 
harmonic. This together with the RMS values of the transverse force 
gives sufficient information about the character! st 1 cs of the transverse 




(M) 



force 
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IV. 01 SCUSS I ON OF RESULTS 



A. BLOCKAGE AND LENGTH-TO-DI AMETER EFFECTS 

Attempts to achieve as high Reynolds numbers as possible in conduct- 
ing wind-tunnel or water-tunnel experiments Invariably give rise to 
wa I l-l nterference effects which, of course, influence whatever measure- 
menTs are desired. There are several blockage correction formulas for 
steady flows which might be used so that the wa I I- interference effects 
might be minimized. Unfortunately none of these formulas could be used 
In the present study for no one has demonstrated that the blockage 
effects In oscillatory flow are identical to those experienced in steady 
flows. 

The blockage ratio D/w, where D is the d iameter of the cylinder and 
w the width or height of the test section, and the length- or span-to- 
diameter ratio, L/D, for the cylinders used In the present study are 
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Achenbach (10) and In some of the experiments of Page and Falkner (II) 
the blockage ratios were 0.166 and 0.185 respectively. The length- to- 
diameter ratio In Page and Warsap f s (12) experiments was 20.2 or 7.88, 
depending on the diameter of the two cylinders they used, as compared 
to 3.33 in the experiments of Achenbach. Evidently, the formulas used 
for steady flow correction effects cannot be applied to oscillating 
flows and that there is not a unique blockage correction for the entire 
period of the harmonic flow. This is evident from the fact that within 
a given cycle the fluid undergoes varying accelerations and velocities 
as the wake width, momentum deficiency, and the wake pressure change 
accordingly. Thus, a blockage correction made for the Instant of maxi- 
mum velocity Is not applicable to the instant at which the maximum 
acceleration occurs. In view of the fact that there are no previous 
investigations, a series of experiments had to be conducted to determine 
the role of blockage in the flow under consideration. For this puroose 
a differential pressure transducer was connected to two pressure taps 
on the same side of the tunnel wall. One of the taps was placed on the 
wall directly above the axis of the test cylinder. The other tap was 
placed 76 cm. (29.92 Inches) to one side of the first tap along a line 
parallel to the flow. A series of experiments were carried out with the 
16.447 cm. (6.475 Inches) cylinder and the differential pressure was 
recorded and compared with the differential pressure obtained from the 
acceleration transducer. Furthermore, to simplify the comparison both 
transducers were calibrated so as to render exactly the same output under 
Identical calibration loads. The results have shown that the two 
differential pressures were nearly identical and that they were certainly 
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within three percent of each other. Often the two traces of two trans- 
ducers were indistinguishably coincident. This somewhat surprising re- 
sult is a clear indication of the fact that the blockage effect in 
harmonic flows is negligible at least for D/w ratios less tnan about 
0.20. Although no special attempt was made to interpret the lack of 
blockage effect in such flows it is believed that the presence of vortices 
on both sides of the cylinder together with the high periods of accelera- 
tion and velocity render the flow relatively more uniform at short dis- 
tances away from the cylinder in the test section. Therefore, for the 
reasons cited above no blockage-effect corrections were appiied to the 
data presented here. It might be of interest to note that had the flow 
been assumed uniform and had the maximum velocity for the largest cylinder 
and the Reynolds number were uc?d to calculate a blockage effect correc- 
tion through the use of one of the existing formulas, one would have 
found that such a correction would have amounted to about six percent. 



B. TRANSVERSE FORCE OATA 
I • Smooth Cy I inde r 

The transverse force coefficients for smooth cylinders are 
plotted as a function of K in figures 8 thrcugh 14. These data 
are also tabulated and presented in Appendix B. It Is evident from tnese 
figures that CLMAX shows some scatter which is in fact smaller than one 
wo* have antici ^ ated in view of the complex texture of the transverse 
force. It is also evident that CLMAX begins to rise rapidly from K » 6 
to 12 and then decreases as K increases. As notec earlier the Reynolds 
number along each plot varies with K and 0. One can designate special 
values of the Reynolds number along tnis curve simply by writing Re =K0 
or by finding the K values which correspond to a particular value of Re 
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Fig. 8 C f (cr.ax) versus K for S = 497. 



39 



SI, ' 



l 1 \ 



I 






V 



i. ! 







Fig. 9 C^(max) versus K for B => 784. 
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Fig. 10 C^(max) versus K for 6 3 1107. 
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Fig. II C^(max) versus K for 6 = 1985. 
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Fig. 12 C^(max) versus K for 0 = 3123. 
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Fig. 13 C^(max) versus K for 8 
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Fig. 14 C^(max) versus K for 0 = 5260. 
45 
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by writing K - Re/8* Prior to doing so, mean lines have been drawn through 
the data as meticulously as possible for each cylinder. Then round numbers 
such as 30,000, 40,000, etc. were chosen for Re and the correspond i ng K 
values were calculated and indicated on each plot. This simple procedure 
resulted In figure 15 in which CLMAX Is plotted as a function of K for all 
smooth cylinders and the identical Reynolds number points are suitably 
connected. It is evident from figure 15 that there is a remarkable and 
heretofore unknown correlation between the force coefficient, Keulegan- 
Carpenter number, and the Reynolds number. The smoothness of the constant 
Reynolds number lines Is another Indication of the consistency of the 
data from one cylinder to another. 

It Is also evicent from figure 15 that for a given Reynolds number 
CLMAX increases at first with K and then decreases gradually. The initial 
rise of CLMAX decreases with increasing Reynolds numbers. In fact for 
Re = 100,000, CLMAX decreases for K values larger than approximately 15. 

For a constant KjCLMAX always decreases with increasing Reynolds numbers. 
The most dramatic change in CLMAX relative to eithe K or Re occurs be- 
tween K equal 8 and 30. The occurrence of the peak CLMAX values in the 
region of K values beTween K equal 10 and 25 Is directly related to the 
shedding of the vortices as will be discussed later. 

The resuits shown in figure 15 are replotted in figure 16 as a 
function of the Reynolds number for constant values of Keu legan-Carpenter 
numbers. This particular plot shows that CLMAX varies significantly 
with Reynolds numbers for Re > 20,000. This observation Is in conformity 
with those made earlier by Sarpkaya and Tuter rsj. For Reynolds numbers 
from about 20,000 to 100,000 the lift coefficient decreases rapidly. 

For larger values of Re, lift coefficient decreases slowly once again and 
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Fig. 15 C. (max) versus K for various values of K and 8. 
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approaches a valje of about 0.3 at Reynolds numbers of about 700,000. 
Evidently, for design purposes the transverse force has significance 
not only because of its magnitude, but also because of its oscillating 
character. Even through the ratio of the transverse force to the in- 
line force is not shown here, a careful examination of the data has 
shown that the transverse forces are of a significant fraction of the 
in-line forces and must necessarily be included in the force calcula- 
tions through a vectorial sum of the two forces. 

As noted above, the oscillatory nature of the transverse force 
may give rise to in-line or transverse hydro-e la*+ic oscillations of 
cylindrical ocean structures whenevpr ine natural frequency of the 
structure falls within approximately +_ 15 percent of the vortex shedoing 
frequency. The in-line oscillations may be excited at natural fre- 
quencies from 2 to 5 times the vortex shedding frequency. It is, there- 
fore, essential that one has definite information about the frequency 
of the transverse force and its dependence, if any, on K and Re. To 
this end the frequency of the first harmonic of the transverse force was 
determined from the analog records for each 0 and the freauency ratio 
defined as 




was calculated. Here f v represents the frequency cf the first harmonic 
and l/T the frequency of the fluid oscillations in the tunnel. The fre- 
quency ratio f^ is shown In figure 17 in such a manner that each curve 
represents the highest values of K and Re below which only frequencies 
indicated can occur. For example, the curve indicated by A shows that 
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Fig. 17 Relative frequency of vortex shedding as a function 
of the Reynolds and Keulegan-Carpenter numbers. 
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to the 1 2 ft of this curve only the frequency ratios of 2 to 8 can occur. 
To the right of curve A, however, frequency ratios larger than 8 are 
possible. The curves for the frequency ratios 3, 5, 7, etc. have not 
been shown in order to simplify the presentat ion. 

It is apparent from the foregoing discussion of the lift coeffi- 
cient, the frequency of the alternating force, and the tabulated data 
that the lift coefficient reaches large values for relatively small 
values of K and Re. It is further apparent that In the region of K 
values (e.g. 8 < K < 15) the relative frequency of vortex shedding is 
not given by an integer number and that there is frequent fractional 
vortex shedding. In other words, the Integer part of f shows the number 
of vortices actually shed while the fractional part indicates that a 
vortex was born, but was never shed. Such a fractional shedding leads 
to f values like 2.1 or 3.7 etc. It also shows that vortices of con- 
siderable strength can grow on one or the other side of the cylinder 
without ever being shed and, thereby, considerably increasing the lift 
force. It is, of course, true that the fractional shedding takes place 
at all values of K and Re; however, the larger the number of vortices 
the smaller is the relative significance of an unshed vortex. It is, 
therefore, apparent that when only two or three vortices are shed, as 
Is the case for K values in the neighborhood of K ^ 15, then the vortices 
which remain attached to the cylinder give rise to significant lift 
forces. It should also be noted, in passing, that the in-line force 
Is signi f icantly effected by the fractional shedding of vortices. 

This becomes particularly evident from the observation of the force 
traces as well as f*om the calculation of the difference between the 
measured force and that predicted from the Mori son equation through the 
use of the Fourier average coefficients. 
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2. Sand-Roughened Cylinders 

The tabulated data are presenrea in Appendix C. The same data 
are also plotted for each cylinder in figures 18 through 25 for two or 
three values of B. 

A quick examination of these plots reveals that 

a) the lift coefficient shows the same + pe of variation with 
K as in the case of the smooth cylinders; 

b) the lift coefficient for rough cylinders is larger than* that 
for the smooth cylinders at the correspond i na Reynolds numbers; and that 

c) the lift coefficient for rough cylinders does not vary with 
B to any appreciable extent. 

A closer examination of the rough cylinder data shows that the effect of 
Re is almost negligible particularly for K values larger than about 15 
within tie range of relative roughnesses and Reynolds numbers tested. 

The transverse force coefficient inevitably exhibits a ‘arger 
scatter than that for the in-line force coefficients because of the some- 
what random nature of the shedding of the vortices. Consequently, it is 
not too uncommon to obtain a variation of 20-2555 for a given K value. 

This fact is of importance in discussing the effect of the Reynolds 
number on the lift coefficient. 

Mean lines drawn through the data shown in figures 22 to 25 for 
each B are combined in figure 26. It is evident from this figure that 
the variation of the lift coefficient from one B to another for a given 
value of K Is smaller than the scatter in any one of the individual 
plots. This shows that within the range of the parameters investigated 
the lift coefficient for rough cylinders does not depend on the Reynolds 
number. Also shown in figure 26 Is the lift coefficient for smooth 
cylinders for B in the range 1000 to 2000. It Is rather surprising 
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Figure 18. C^(max) versus K for k/D = 1/370 
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Figure 19. C L <max) versus K for k/D = i/370 
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Figure 20. C^(max) versus K for k/D = 1/370 
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Figure 21. C^(max) versus K for k/D = 1/370. 
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that the smooth cylinder data at relatively low Q values form more or 
less the upper limit of the rough cylinder data. Thus, the following 
important and somewhat surprising conclusion has been reached: for 
smooth cylinders, the lift coefficient depends primarily on Re and very 
little with K for Re larger than about 20,000, For rough cylinders, on 
the other hand, the lift coefficient ceases to depend on Re and varies 
primarily with the Keu I egan-Carpenter number. Finally, the lift co- 
efficients for rough cylinders at high Reynolds numbers ar£ almost 
identical with those for smooth cylinders at very low Reynolds numbers. 
Even though the reasons for the foregoing are not quite clear, the 
consequences of this experimental finding for model testing purposes 
are rather obvious. Suffice it to say that roughness has a profound 
effect on the lift coefficient. This is partly because of the fact that 
it gives rise to vortices of larger strength and hence to larger lift 
forces, but also because the lift force for a rough cylinder is a larger 
fraction of the in-line force. For. example, at a Reynolds number of 
700,000 the drag coefficient for a roi . h cylinder with a relative 
roughness of 1/370 is about 1,4 for K = 50 and the lift coefficient at 
the same correspond ing K and Re values is about 1.0. Thus, in all 
designs the vectorial sum of the in-line and transverse force will have 
to be taken into account. Otherwise the design force will be grossly 
underestimated. 

As noted earlier, the alternating nature of the transverse force 
is as important as its magnitude. It is for this reason that the fre- 
quency of the alternating force has also been calculated, as in the 
case of smooth cylinders, and tabulated (see Appendices B and C). A 
close examination of the frequency ratios shows that f r /K remains 
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essentially constant at a value of about 0.22. To be sure, there are 
variations from one cylinder to another and from a given combination 
of Re and K to another one. Neverthe less, the Strouhal number is fairly 
constant for all roughnesses, relative amplitudes, and Reynolds numbers 
larger than about 20,000. 

In the foregoi ng, the major features of the transverse force act- 
ing on smooth and rough circular cylinders in harmonic flow have been 
described. It has been shown that the lift force is a significant 
fraction of the total force acting on the cylinder and that it is in- 
tensified by roughness. It appears that the larger shear stresses on 
rough cylinders give rise to vortices of stronger circulation and 
thereby increase the lift force. It would, therefore, be extremely 
interesting to measure the pressure distribution around a rigidly 
mounted rough cylinder and compare the results with those obtained with 
a smooth cy I i nder. 

Finally, a few words are necessary regarding the application 
of the results presented here to the structures in tne ocean environ- 
ment. It is realized at the onset that the flow induced by waves is 
not two-dimensional and that the velocity decreases with depth and the 
fluid particles move on elliptical paths. The variation of the velocity 
along a vertical pile invariably results in a decreased correlation in 
the vortex shed. This in turn results in a reduced transverse force. 
Thus, the results presented herein must be applied with the full reali- 
zation of the foregoing facts and must be regarded as the maximum upper 
limit of the transverse force that could be generated in wavy flows. 
Evidently, at great depths where the change of velocity is very small 
along a vertical pile, relatively large lift forces must be expected. 
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At elevations closer to the ocean surface, however, the rapid changes in 
velocity may decrease the correlation length and thus give rise to 
smaller transverse forces. It has been assumed that these structures 
have been rigid and are rigidly mounted and that there is no motion due 
to the forces acting on them. In general most members of such a structure 
are elastically mounted and that they may undergo self-excited oscilla- 
tions under the action of waves and/or currents. The foregoing discus- 
sion points out not only the regions c* application of the results pre- 
sented herein, but also the many unsolved complex problems which are 
encountered in the real ocean environment. It is, therefore, suggested 
that further investigations be undertaken to study these hydro-elastic 
oscillations of smooth and rough cylinders in harmonic flows, in-line 
and transverse forces on yawed cylinders, and finally the effect of 
the wave-current combination on rigid or elastic, smooth or rough 
cylinders, near to or away from a wall, and normal or inclined to a two- 
dimensional harmonic flow. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The extensive investigation of the transverse forces on smooth and 
rough circular cylinders in harmonic flow warrants the following conclu- 
sions: 

1. The transverse force coefficient depends in general on both the 
Reynolds and Keu I egan-Carpenter numbers. The variation of CLMAX with K 
and Re depends not only on the particular value of Re but also on the 
relative roughness k/D; 

2. For smooth cylinders, CLMAX does not vary appreciably with Re 
for Re smaller than about 20,000. In this range the variation of CLMAX 
with K is such that It reaches a maximum at about K - 15 and then de- 
creases as K Increases. For Reynolds numbers between 20,000 and 150,000, 
CLMAX varies rapidly with Re and then gradually reaches a value of about 
0.25; 

3. The relative frequency of vortex shedding, when expressed in 
terms of the Strouhal number, shows that for Reynolds numbers smaller 
than about 55,000 the Strouhal number remains essentially unchanged; 

4. For rough cylinders, the lift coefficient has several heretofore 
unknown and practically significant characterist ics. The results show, 
within the range of the parameters tes'^d, that CLMAX does not depend 

on Re. The resulting CLMAX distribution is, surprisingly enough, very 
close to that obtained with the smooth cylinders at very 'ow Reynolds 
numbers. 

5. The Strouhal number for rough cylinders remains nearly constant 
for all Reynolds numbers at about 0.22 with the possible exception of 
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those at very low Reynolds numbers. It should be emphasized that one 
could not have reached this conclusion had one not carried out the 
painstakingly difficult series of experiments described in this thesis 
6. It Is recommended that pressures be measured around smooth and 
rough cylinders and that the effect of the self-excited oscillations 
of cylinders on the In-line an* transverse forces be examined. It is 
hoped that through experiments of Intrinsic quality one will be able 
to isolate the Influence of the major factors involved and arrive at 
sound rules for design of structures in the ocean environment. 
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SPECIAL MOTES FOR APPENDICES B AND C 



The data given in the following tables are identified with the 
diameter of the cylinder, relative roughness, and the frequency 
parameter . 

The following symbols have been used: 



DIAM - diameter of the test cylinder 

k - height of sand roughness 

D - diameter of the test cylinder 

2 

6 - frequency parameter (D /vT) 

REYN - Reynolds number (U D/v) 

1 m 

UMTOD - Keu I egan-Carpenter number (U^T/D) 
CLMAX - maximum lift coefficient (see eqn. I) 
FR - frequency ratio (see eqn. 3) 

3T - c *rouhal number (see eqn. 4) 
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APPENDIX - 8 TRANSVERSE FORCE DATA FOR SMOOTH CYLINDERS 
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APPEND I X-C TRANSVERSE FORCE DATA FOR ROUGH CYLINDERS 
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